Systematics of geometric scaling 
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Using all available data on the deep-inelastic cross-sections at HERA at a; < 10~^, we look for 
geometric scaling of the form ''(r) where the scaling variable r behaves alternatively like log Q^— 
XY, as in the original definition, or log Q^—X\/Y, which is suggested by the asymptotic properties of 
the Balitsky-Kovchegov (BK) equation with running QCD coupling constant. A "Quality Factor" 
(QF) is defined, quantifying the phenomenological validity of the scaling and the uncertainty on the 
intercept A. Both choices have a good QF, showing that the second choice is as valid as the first 
one, predicted for fixed coupling constant. A comparison between the QCD asymptotic predictions 
and data is made and the QF analysis shows that the agreement can be reached, provided going 
beyond leading logarithmic accuracy for the BK equation. 



1. Geometric scaling Q is a remarkable empirical 
property verified by the data on the high-energy deep- 
inelastic scattering (DIS) cross-sections a'^ p. It has been 
realized that one can represent with reasonable accuracy 
the cross-section by the formula 
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where Q is the virtuality of the photon, Y the total rapid- 
ity in the 7*-proton system and oc e^^ an increasing 
function of Y . The value found for A ~ 0.3 has been con- 
firmed _by the well-known Golec-Biernat and Wiisthoff 
model |3| where geometric scaling has been explicitly as- 
sumed in the parametrisation. 

The property also observed in DIS on nucleus (3| 
and diffractive processes has been intimately related 
[3 to the concept of saturation 0], i.e. the behaviour 
of perturbative QCD amplitudes when the density of 
partons becomes high enough to exercise the unitarity 
bound. Indeed, there has been many theoretical argu- 
ments to infer that in a Y, domain where saturation 
effects set in, geometric scaling is expected to occur. 
Within this framework, the function QsiY) can be called 
the saturation scale, since it delineates the approximate 
upper bound of the saturation domain. 

Following a theoretical approach, it has been possible 
to derive the property from the nonlinear Balitsky- 
Kovchegov (BK) equation which represents the "mean- 
field" approximation of the evolution equation for high 
energy (high density) QCD amplitudes. This equation 
is supposed to capture essential features of saturation 
effects. In the 1-dimensional approximation it reads for 
the density of gluons with transverse momenta k in some 
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dyUiL, Y) = axi-dL) Af{L, Y) - aAf^L, Y), (2) 



where L — \og{k^ /k"^), with fcp being an arbitrary con- 
stant. In Eq.(|21), the coupling constant will be con- 
sidered alternatively in the following as fixed (we shall 
use d = 0.15), or running such that d{L) = 1/bL ,6 = 
{llNc-2Nf)/l2Nc. The kernel wiU also be either taken 
with leading logarithm (LL) accuracy as xil) = 
2?/'(l) — -0(7) — -0(1 — 7) , or at next-to-leading logarithm 
(NLL) accuracy (see e.g. |l[lS[I2). 

The key ingredient to theoretically prove geometric 
scaling of the asymptotic solutions of the nonlinear equa- 
tion ((21) is the travelling wave method [l^ . Indeed, the 
BK equation admits solutions in the form of travelling 
waves Af{L — VgdY). L has the interpretation of a space 
variable while t — dY, interpreted as time, is an increas- 
ing function of rapidity Y. Vg is the critical velocity of 
the wave, defined in this case as the minimum of the 
phase velocity. The travelling-wave solution for the quan- 
tity Af can be easily translated to the property (|2Jl (at 
least assuming negligible quark masses) since it yields 
N [k"^ / Ql{Y)) . Hence the asymptotic travelling- wave so- 
lutions of the BK equation satisfy geometric scaling (up 
to subdominant scaling violations which we will not anal- 
yse in the present work). Similar results can be obtained 
|l3l Il4l | from an approximation of the linear equation 
with absorptive boundary conditions. The main theoret- 
ical results for the asymptotic saturation scale are the 
following: 

logQl = d^{Y - Fo) - log(r - Fo) 
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(a) Quality Factor as a function of A. 



(b) Scaling curve. 



FIG. 1: Geometric Scaling in Y 
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for running coupling. 7c is the solution of the implicit 
equation xilc) = JcX'ilc) while = —2.338 is the first 
zero of the Airy function. Yq is an arbitrary constant 
which may parameterize the unknown "non-universal" 
preasymptotic contributions, depending on the initial 
conditions. 

The questions we want to ask deal with the confronta- 
tion between the empirical formula |(2J) with the asymp- 
totic predictions for the saturation scale (|3I4|1 . They may 
be expressed as follows 

• Is geometric scaling, demonstrated for log(Qg) ~ Y 
in agreement with the leading term in remains 
vahd for \og{Q'j.) ~ Y^^^ as suggested by (gj? 

• Are the theoretical factors and subasymptotic 
terms given in Eas. (|3l4|l visible in the data? 

• Are NLL effects in the theoretical predictions phe- 
nomenologically important? 

• What is the role of "non-universal" terms? 

2. We focus our analysis on the measurements of the 
7*p cross-section a'' p — (47r^Q;e/Q^)^2 for which geo- 
metric scaling is predicted. The latter being valid at 
high-energy, we shall restrict ourselves to a: < 0.01 for all 
values of Q^. Within that range, we shall use all avail- 
able data (404 points from E665 0, HI 0, NMC [13 
and ZEUS [H). 



In order to get a quantitative answer, we shall intro- 
duce an "estimator" or quality factor (QF) for determin- 
ing the scaling quality, whose working definition is the 
following. Given a set of points {xi,yi,fi = f{xi,yi)) 
and a parametric scaling law t{x, y; A), we want to deter- 
mine if for some values of A, /(x, y) can be considered a 
function of t only. To achieve this we define a QF which 
is large when the points (r^ = T{xi,yi]X), fi) "lies on a 
unique curve". To quantify this still ill-defined concept, 
we shall first rescale the set (t^, fi) into (wi, Vi) such that 
Q < Ui,Vi < 1 and assume that the u^'s are ordered. We 
then introduce 
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This definition for the quality factor obviously achieves 
what we want: when two successive points are close in 
u and far in v, we expect them "not to lie on the same 
curve" and, indeed, they give a large contribution to the 
sum in |(SJ), leading to a small quality factor. The con- 
stant £^ in ((Sj is a small number (we have taken e = 1/n 
with n being the number of data points) which prevents 
the sum from becoming infinite when two points have the 
same value for u. Finally, by studying the dependence on 
A of Q (as we shall see, it usually shows a Gaussian peak) 
we can determine the best choice for the parameter in the 
scaling function r and its uncertainty. 

3. In order to answer the first question concerning the 
compared validity of geometric scaling in Y vs. Vy, we 
performed the QF analysis in both cases. In FigQ] we 
show the data and QF for the usual geometric scaling 
definition with Y. In Fig^(b), is displayed the data plot 
with scale redefinition. The quality of scaling on this plot 
can be read from FigQKa), where the scatter plot of the 
QF is displayed together with a Gaussian fit of the peak. 
As obvious from Fign}(b), the peak value larger than 0.1 



3 




showss a good scaling and a Gaussian fit of the bump 
gives the best value and error for the lambda parameter 
A = 0.321 ± 0.056 . Note that the results displayed on 
Figncan be considered as an update for the usual geo- 
metric scaling tests and a quantitative determination of 
the scaling parameter and its attached uncertainty. The 
new point is that geometric scaling is also verified, and 
with the same level of quality than previously, when the 
saturation scale is chosen as a function of ^Jy . This is 
manifest on FigEl where comparable QF heights at the 
peak and width attest of the quality of scaling. The peak 
value getting larger than 0.11 givess a good scaling. The 
Gaussian fit of the bump gives A = 1.621 ± 0.254 . 
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FIG. 3: Quality Factor for tlie LL prediction. 



4. Having verified geometric scaling, as suggested by 
the asymptotic solutions of the QCD equations with both 
fixed and running coupling constant, it is tempting to go 
one step further and to ask whether the theoretical for- 
mulae (|3I4|I with their theoretically predicted parametri- 
sations are valid? 



For that sake, we have considered first the LL pre- 
diction with fixed coupling Q with the first subleading 
correction, i.e. of the form log((5^) oc AF — /j,log(F) (see 
Fig. 121). The predicted values from are represented 
by a point in the A, /i plane. It happens that a QF of 
'-^ 0.055 can be attributed to this parametrisation, as 
shown on the same figure by the equi-QF curves. In this 
case changing Yq or adding the third term of Q does not 
improve the QF's. Hence the asymptotic LL theoretical 
prediction remains only marginally verified at present en- 
ergies. 

Let us now turn to the analysis of the 3 different the- 
oretical NLL schemes that we consider. As now well 
known, NLL corrections to the LL kernel ^] have to be 
embedded in a resummation scheme in order to cancel 
spurious singularities. We consider asym pto tic predic- 
tions igj for the so-called S3, S4 scheme s |loll and CCS 
scheme 11] which were recently derived [l9|. 

On Fig^ we show the resulting analysis. The Qual- 
ity Factor for the 3 schemes is presented as a function 
of Iqi which parameterizes in the non-universal con- 
tributions. There always exists a range of Yq for which 
the QF goes beyond 0.1 and thus leads to an accept- 
able scaling. This has to be contrasted with the situa- 
tion with LL kernel. However, if Yq is too negative, the 
relevance of geometrical scaling in \/Y is questionable. 
Indeed t = ^Y -Yq - + Y/2y^\ and thus the 

scaling variable is Y with a non universal A = l/2-y/|lo| 
parameter. For that reason the S4 scheme seems favoured 
by the analysis and the figure displays the resulting scal- 
ing curve for Yq = —5.5. 

5. Using the preceding results based on the estimate of 
Quality Factors for scaling properties, it is now possible 
to answer the questions asked in the introduction. 

• Geometric scaling is definitely as valid for the 
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(a) Quality Factor as a function of Yq 
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(b) Scaling curve for the S4 scheme. 



FIG. 4: Geometric Scaling at NLL accuracy 



choice log((5g) ~ y^/^ as it is for the original 
log(Qs) y suggestion 

• The theoretical predictions JSJ based on the BK 
equation at leading log accuracy, i.e. with fixed 
coupling and LL kernel are only marginally verified 
with a QF of ~ 0.055. 

• At next-to-leading accuracy, the theoretical predic- 
tions l@J on scaling may be satisfied. There exists 
preferable NLL schemes e.g. the S4 scheme [foj . 

• Parametrising "non-universal" terms, i.e. depend- 
ing on initial conditions by a constant Yq in the 
rapidity evolution, we find that they do not play 
an essential role at LL level, while they contribute 
to the quality of the scaling at NLL level. 

The generality of the Quality Factor method indicates 
its interest for a quantitative evaluation of good scaling 
properties beyond the specific applications we focussed 
on in the present work. It can be used to check empirical 
scaling-law proposals as well as to quantitatively evaluate 
the "distance" between given theoretical predictions with 
data. In the present application, it leads to the conclu- 
sion that geometrical scaling is a well verified empirical 
property of deep-inelastic cross-sections on the proton. It 
also shows that the theoretical predictions based on the 
QCD saturation mechanism are consistent with the ob- 
servation, but requires to go beyond leading logarithmic 
accuracy and the asymptotic terms. 

Finally, let us quote that recent high-energy predic- 
tions 20] of diffusive scaling, namely cr^ ^/VY — o'(t) 
with T = \q%[Q'^ I Ql) I \/Y , have a QF of 0.05, indicating 
that it might require higher energies. 
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